Early out-of-equilibrium beam-plasma evolution 
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We solve analytically the out-of-equilibrium initial stage that follows the injection of a radially 
finite electron beam into a plasma at rest and test it against particle- in- cell simulations. For initial 
large beam edge gradients and not too large beam radius, compared to the electron skin depth, the 
electron beam is shown to evolve into a ring structure. For low enough transverse temperatures, 
the filament at ion instability eventually proceeds and saturates when transverse isotropy is reached. 
The analysis accounts for the variety of very recent experimental beam transverse observations. 

PACS numbers: 52.35.Qz, 52.40.Mj, 52.65.Rr, 52.57.Kk 



Beam-plasma interactions have recently received some 
considerable renewed interest especially in the relatively 
unexplored regimes of high beam intensities and high 
plasma densities. One particular motivation lies in the 
fast imition schemes (FIS) for inertial confinement fu- 
sion These should involve in their final stage the in- 
teraction of an ignition beam composed of MeV electrons 
laser generated at the critical density surface with a dense 
plasma target. The exploration of the electron beam 
transport into the overdense plasma is essential to assess 
the efficiency of the beam energy deposit. In this mat- 
ter, transverse beam-plasma instabilities could be partic- 
ularly deleterious in preventing conditions for burn to be 
met. Experimental observations recently undertaken in 
conditions relevant to the FIS have either shown some 
transverse microscopic filamentation of electron beams 
[2] or some transverse, predominantly macroscopic, beam 
evolution into a rin^ structure 0, Q or a superposition of 
those effects ^ 1^ ? with filaments standing out from 
a ring structure, in a scenario similar to Taguchi et al.^s 
numerical simulations 7]. Weibel instability 8] is com- 
monly invoked to account for these phenomena, but it 
is sometimes difficult to find any clear univocal evidence 
supporting this. Moreover the fact is that, whereas most 
theoretical and some computational studies are devoted 
to the linear regime of instabilities originating from cur- 
rent and charge neutralized equilibria, the physics of the 
fast ignition is intrinsically out-of-equilibrium. 

In this Letter, we shall consider the out-of-equilibrium 
initial value dynamical problem taking place when a radi- 
ally inhomogeneous electron forward current is launched 
into a plasma and is still not current compensated. We 
shall focus on this early stage where collisions may be ne- 
glected. Ions will be assumed to form a fixed neutralizing 
background. In order to simplify both the analysis and 
the numerical PIC computations, we shall consider the 



system to be infinite along the beam direction z. We shall 
remove any z dependance by assuming also that plasma 
density Upe is uniform and constant. At time t = 0, a 
relativistic electron beam is switched on in the plasma. 

Maxwell equations are linear and can thus be solved for 
all time to give the electromagnetic fields as functions of 
the sources, namely beam and plasma current densities, 
jb and ipe. We get rotB = /io(jpe + jb) + l/c^dE/dt 
and rotE = —d'B/dt. The electron plasma current jpe = 
—eUpeVpe is initially vanishing and may be approximated 
by linear fluid theory in the initial stage yielding 
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with ujpe = ^J Tipee^ I vfieS^ the plasma pulsation. We 
Fourier decompose any field g through g{r,0,t) = 
Xlm 9^^\'^'> t) exp {imO) and proceed to a Laplace trans- 
form in time g'^^^ (r, s) = ^-stg{m) (^^^ Elimi- 
nating the electric field components. Maxwell equations 
in cylindrical geometry yield inhomogeneous wave equa- 
tions with sources for the magnetic field components. In- 
troducing the operator jCn[fA such that 
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defining a = y 5^ + ^pe/^ neglecting the initial vaL 
ues of the e.m. fields, the wave equations read 

^(0) _ „ dfbJ 



for m = 0, £1 [ia] 
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form 7^ {),iCm[i(j] (rB^"^^^ 
with, for any m. 
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Let us make the fohowing general statements: Because 
they are linear, Maxwell equations do not enable spec- 
tral changes. If the beam is sufficiently weak, so that 
the fluid approximation for the bulk plasma remains ap- 
proximately valid, mode transfers will originate from the 
beam particles equations of motion. Consequently, if the 
initial beam is rigorously both rotationally and axially in- 
variant (on m = 0), it will remain so for all times. Then 
there are no sources to feed the triad {Br, Bz, Eq) that 
remains vanishingly small. This is an invitation to focus 
first on the m = evolution. 

We consider initial beam density and velocity of 
nhQ>{r) = nbon{r) and VbOz{r) = VbOzv{r), with f = r/n. 
Let us introduce here the beam radius 9], the elec- 

their ratio r] = r^/Xs, and 



tron skin depth Xg 
let us define a = 

tial relativistic Lorentz factor 70 = (l 

-1/2 



nbo/ripe, Po = Vboz/c and the ini- 
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(1 — [3nv'^) ^'^ — 7o(^). The Green function g(r \ a) 
jinj solving Ci [ia] g = —S {r — a) is readily computed as 
g{r I a) — h {ar'^) Ki (ar^) with = min(r, a) and 
= max(r, a). The general solution of Eq. is then 



B^^\r,s) = -fio I ah{ar<)K^{ar>)^da. (6) 
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Let us consider the response Bq^^ to the initial beam 

current = —en\)QV\)QzjH(f) that is switched on at 
time 0. Here H denotes the Heaviside step function and 
j(r) = n{r)v(r). This gives J^q^^^^.s) = -enboVbOzj{r)/s. 
Eq. admits then a solution in separate variables. This 
makes Laplace inversion easier giving 



—B^^-^ (f, t) = af3ou;per]F{r, rf) [1 - cos {ujpet)] 
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where the radial information is contained into 

F{r^r]) = Ki{r]r) / uIi{r]u)j\u)du-\-Ii{r]r) 
Jo 



We now wish to compute the beam evolution under 
the previous self- fields {TJ, which was ignored in previ- 
ous studies Let us assume that the beam can be 
treated as a cold fluid and write the m = fluid equa- 
tions, dropping the m = superscripts. 
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Let us write Uh = ^6o+^6i(^, t)^ v^r = '^6ir(^, t) and v^z = 
VbOz + Vbiz{r^ t). Considering a as a small parameter, we 
get a natural hierarchy: first order terms should be of 
order a, second order terms of order and so on. Let 
us explicit first order fluid equations. The radial electric 
contribution being negligible (|Tn|) . Eq. (|T^ gives 
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which yields, using vi}ir{r, t = 0) = 0, 
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The first order conservation equation is 
dribi 
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Then, integrating dE^^ /dr 
(fTJ immediately gives 



dB^Q^ /dt and using Eq. 



uKi{r]u)j\u)du. 
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Using (|T^ . this gives, with nbi{r^t = 0) = 0, 

^bl{r,t) ^2- A 2 .2 , / .^ 
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where ni is the radial function 
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Finally, the m = radial electric field component satisfies 
the wave equation duEr^^ ^u'^^Er^^ = —l/sodtj^^- It is 
easy to check that its initial behavior is given by 



Eq. (|T7|) shows that ribi has a secular behavior. Thus the 
present analysis breaks when ribi/ribo = 0{1)^ namely 
roughly for copet ~ In order to put this more 

precisely, we shall study the radial behavior ni. Let us 
consider some initially monokinetic beam ('^ = 1) having 
density functions of the form 



n(f) = (tanh [2/^(1 - f)] + 1) /2. 



(19) 



— E^^\r, t) = aul^n{f) I cos [upe{t - r)] %|,(r, r)(ir. 
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(10) 

For the problem under consideration, the beam to plasma 
density ratio a is typically a small parameter. Eq. ([Tn|) 
will be second order in a. 



This enables the study of the influence of the beam edge 
gradients as z^: = —n'{f = 1). Fig. Q displays various 
n profiles and their associated first order perturbations 
ni. This shows that, within a given initial lapse of time, 
the natural evolution of the system tends to increase the 
beam density around some radius below r^. This favors 
the formation of a ring structure for the beam density. 
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n(r) 




FIG. 1: Functions n(f) (tanh [2k(1 - r)] + 1) /2 and ni(f) 
for ?7 = 10 and k = 10 (bold), = 5 (dashed line) and k = 2 
(plain line) for a monokinetic beam. 

that is all the sharper and all the closer to that the 
initial radial beam gradients are high. The emergence of 
this beam ring formation can be already inferred from 
the time evolution of test electrons within the azimuthal 
magnetic self-field corresponding to some given initial 
profile n of the form (|T^ as shown in Fig. [3 The caus- 
tics pattern signals there a cusp formation in the radial 
beam density. It is important to note that, due to Eq. 

this behavior is very dependent on the initial beam 
profile. Indeed, we observed no such caustics pattern nor 
any visible evolution towards a ring structure for smooth 
Gaussian initial beam profile, but rather radial focusing 
prior to the filament ation onset. 
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FIG. 2: Radial trajectories of test electrons for an initial beam 
profile given by Eq. (|19|) with 77 = At = 10. 

In order to assess the validity of the first-order ana- 
lytical results presented above and get an insight into 



the longer time evolution of the beam-plasma system, 
we performed particle-in-cell (PIC) simulations using the 
code C ALDER 12] in 2-1/2 dimensions (x^y^Vx^Vy^Vz). 
All species, namely plasma ions and electrons and beam 
electrons, are described as particles. Beam electrons are 
injected at t = in a plasma without current compensa- 
tion and As = 0.05 /im. Fig. 01 presents the early time 
evolution of the beam radial and poloidal mean veloci- 
ties. The initial beam density profile was given by Eq. 
(|T9j) with ?7 = 10 and Hi = 2. The beam was monoki- 
netic with 7o = 15 and a was equal to 0.03. For beam 
radial velocity, this figure shows a nice agreement with 
the analytical result (P3|) . As for beam poloidal veloc- 
ity, it is initially vanishing and its m = component 
does remain so. However, poloidal symmetry eventu- 
ally breaks due to the arising of filamentation instabil- 
ity. This takes place after a short transient during which 
plasma backcurrent grows. Then the average poloidal 
velocity start to grow exponentially, with a growth rate 
that nicely fits the linear filamentation instability one, 
given by 7/ = Poy^a/jocVpe [Ulil. 



Numerical 
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FIG. 3: Early time evolution of the modulus of the radial 
average of the transversal beam velocity {Pr = Vr/c and pe = 
vo/c) in lin-log scale for a = 0.03 and 77 = 10, a>: = 2 in (fTo)) . 
Beam and plasma transverse temperatures are IkeV. 

Fig. 21 presents the longer term evolution of the trans- 
verse components of beam velocity for a monokinetic 
beam with a larger value of a {a = 0.15) and smaller 70 
(70 = 3). As previously, there is an initial phase, between 
t=0 and 20uj~J-, during which radial velocity grows fast 
and poloidal velocity remains small. For cjpet o:^ 20, we 
can see on the inset of Fig. ^Ithat beam density presents a 
clear ring structure at its edge. When the beam current 
is partially neutralized, filamentation instability starts, 
breaking the initial azimuthal system symmetry and pro- 
ducing the exponential growth of poloidal beam velocity. 
When the instability saturates {ujpet 35), the mag- 
nitudes of both components of the transverse velocity 
are similar: transverse isotropy is reached. For 77 larger 
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(not shown here), the relative thickness of the initial ring 
diminishes and once filament at ion saturates, the initial 
structure becomes almost undetectable. 
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FIG. 4: Same as Fig. Hfor a = 0.15, 70 = 3, ry = 10 and 
K, = 20. Transverse snapshots of beam density are included. 
Beam and plasma transverse temperatures are 5 keV. 



a laser pulse over a solid target may present divergencies 
as large as 17°. The origin of these large divergencies is 
not clear. We performed simulations with beams having 
an angular divergency of 15° and two values of emittance. 
For zero emittance (laminar beam) a dim ring structure 
appears at short times (Fig. ISI^a)), then filamentation 
takes place (Fig. [S^b)). The large radial velocity pro- 
duces a fast coalescence of the filaments along the radial 
direction, resulting in a star-like density pattern. In the 
high emittance case (Fig. EJc-d)), no ring structure is ap- 
parent. Moreover, the large transverse temperature pre- 
vents the onset of the filamentation instability |lj, |l^ • 

In conclusion, our analysis has shown that, depending 
on its initial radial density and velocity distribution, the 
shape of an electron beam propagating in a plasma may 
evolve into a transient ring structure. This results from 
the natural evolution of the system and not from the 
usually invoked Weibel instability. If its transverse tem- 
perature is low enough, filamentation instability eventu- 
ally proceeds. The observation of the ring structure is 
favored by sharp beam edges and not too large beam 
radius (compared to the electron skin depth). It is not 
generic which may explain the variety of experimental 
observations 0,11,13,1111. 

Discussions with A. Bret are gratefully acknowledged. 
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FIG. 5: Spatial distribution of the beam density at two times 
for beams with angular divergence of 15°. Zero emittance 
(a-b) and large emittance (0.25 \im) cases (c-d) are shown. 

Finally, the evolution changes significantly for beams 
with finite initial angular divergence. Experimentally, it 
was found [l5| that electron beams created by focalizing 
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